function [c,y] = eqZD_func_con(x,i_m,i_D,lambda_func,alpha1,alpha2,alpha3)
%===============================
%equilibrium given D rate and im
%===============================
y(1)=alpha1*lambda_func(x(1))+alpha3*lambda_func(sum(x))-i_m;
y(2)=alpha2*lambda_func(x(2))+alpha3*lambda_func(sum(x))-i_D;
c=[];

end